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SYMBOLS 


c 

F 

f 

G(r, t; r s ,t s ) 

G(f,i;f s ,i s ) 
Gh(?; ?,) 

G2 d(x,y;x s ,y s ) 

Hj(*) 

l 

J n (a:) 

k 

M 

m 

N 

n 

n 

P 

P 

R 

rt 

? s 


speed of sound in flow 

source term in acoustic equation 

source term in Helmholtz equation 

Green’s function for acoustic wave equation in uniform 

subsonic flow 

Green’s function for acoustic wave equation without flow 
Green’s function for subsonic convected Helmholtz equation 
Green’s function for subsonic convected Helmholtz equation 
in two dimensions 
Hankel function, Jo(z) + «Yo(x) 

/=T 

Bessel function of the first kind of order n 
wave number 

wind tunnel Mach number, Uqo /c 
harmonic number of acoustic frequency 
number of blades 

normal to duct surface, positive toward duct center 
harmonic number of loading or thickness function 
acoustic pressure 

coefficient of Fourier component of acoustic pressure 

radial distance 

rotor radius 

position 

source position 
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t s 

U 00 
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x,y,z 
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Y n (x) 
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duct surface 
time 

source time 

velocity of uniform flow in duct 
duct volume 

coordinates in fixed system 

coordinates in moving system 

Bessel function of the second kind of order n 

boundary condition coefficient 

boundary condition coefficient 

Dirac delta function 

three-dimensional Dirac delta function 

altitude angle of observer to rotor plane 

distance along blade chord 

average density in flow field 

azimuth angle of observer 

acoustic frequency 

rotational frequency of rotor 



THE INFLUENCE OF A WIND TUNNEL ON HELICOPTER ROTATIONAL 


NOISE: FORMULATION OF ANALYSIS 
Marianne Mosher 
Ames Research Center 


SUMMARY 


This paper describes an analytical model that can be used to examine 
the effects of wind-tunnel walls on helicopter rotational noise. First, a com- 
plete physical model of an acoustic source in a wind tunnel is described and a 
simplified version is then developed. This simplified model retains the impor- 
tant physical processes involved, yet it is more amenable to analysis. Second, 
the simplified physical model is modeled as a mathematical problem. An in- 
homogeneous partial differential equation with mixed boundary conditions is 
set up and then transformed into an integral equation. Details of generat- 
ing a suitable Green’s function and integral equation are given in appen- 
dixes. Last, the equation is discussed; it is also given for a two-dimensional 
case. 


INTRODUCTION 


It has become important to know the noise levels that a new helicopter 
will produce before it is put into production. Rules of the International Civil 
Aviation Organization (ICAO) and possible future rules of the Federal Aviation 
Administration (FAA) place limits on the allowable noise, measured in terms of 
effective perceived noise level (EPNL), that a helicopter may produce. To meet 
these rules, both EPNL and more detailed noise measurements of a helicopter 
need to be known before manufacturers produce new or modified helicopters. 
The current analytic prediction methods for the noise levels are completely 
inadequate for newly designed helicopters. The totally theoretical methods 
available do not give predictions that correlate well with measurements (for 
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example see the comparisons in refs. 1 and 2). Using measured air loads 
in the formulations improves the comparisons, but does not yield satisfactory 
results either (refs. 3 through 5). Empirical noise prediction methods are 
used to give better correlation with measured noise levels (ref. 6) . However, 
an empirical method will only yield good results for helicopters that are very 
similar to those used to produce the empirical predictions. A newly designed 
helicopter will produce noise different from that of earlier helicopters; moreover, 
a new helicopter will have to be quieter than many existing helicopters in 
order to meet proposed noise rules of the FAA and ICAO and community noise 
standards. 

Measuring the noise levels produced by small-scale and prototype models 
of a new-design helicopter is another way to estimate the noise The noise 
levels of interest are those produced during flyover, takeoff, and landing. 
Measurements made in any wind-tunnel test must be translated to the cor- 
rect reference frame for a flyover measurement Measurements from a full- 
scale prototype must be corrected for the reverberant effects < f the vind-tunnel 
walls, for there are no anechoic wind tunnels large enough io lest a full-size 
helicopter Measurements from a small-scale model must be io. lotted for scale 
effects and may have to be corrected for reverberant wind-tunnel effects as 
well. 


PHYSICAL MODEL OF PROBLEM 


The problem of predicting the acoustic field of aeroacoustu M,.irces in a vind 
tunnel is complicated A model that incorporates all of the ph^ steal phenomena 
is very elaborate. First, the source and its acoustic radiation must be modeled 
in a flow field The drive system of the wind tunnel also intioduces a noise 
source. The presence of the wind-tunnel walls may have an effect on the source 
radiation (or strength). The walls produce a complicated semi-reverberant sound 
field by reflecting the sound. The boundary layer on the wall may alter the 
acoustic reflection. Also, wind-tunnel walls are a thin shell, therefore, an acoustic 
field inside the wind tunnel is coupled to structural modes of the wind tunnel 
which are also coupled to sound radiation from the wind-tunnel walls to the 
outside. Finally, there is the geometry of the wind tunnel to consider The mean 
flow may vary because of changes in the cross-sectional area, thus complicating 
the acoustic field. Obstacles, such as vanes, screens, and drive fans, as well as 
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changes in cross-sectional area, reflect acoustic energy in the wind tunnel. If the 
tunnel is an open circuit configuration, acoustic pressure at the ends of the duct 
couples to the acoustic field outside the tunnel radiating away from the wind 
tunnel. 

Modeling all of these phenomena requires a complicated system of coupled 
partial differential equations with each equation modeling a different aspect of 
the problem. The following describes the complexity of this type of model. 
The Ffowcs- Williams and Hawkings equation with a moving observer models 
the aeroacoustic source and its radiation to a surface just outside the boundary 
layer of the wind-tunnel wall. Acoustic pressure and velocity couple at the 
boundary layer to another equation describing acoustic transmission through 
the boundary layer. At the wall, the acoustic pressure and velocity couple 
to the wall velocity and acceleration. A fourth-order structural equation 
governs the vibration of the walls. Velocity and acceleration on the outside 
of the wall couple to the acoustic field outside. If the wind tunnel is an 
open circuit, the acoustic pressure field at the inlet and exit of the tunnel 
couple to the outside acoustic field. Outside the tunnel, the acoustic field 
can be modeled as a simple potential field with a radiation condition in the 
far field. At present the above model is impossible to solve, except for 
possibly the most elementary case involving very simple sources and tunnel 
geometries. 

This paper presents an investigation of the sound field inside the wind tunnel 
relatively close to the source, so a simpler model, which retains the essential physi- 
cal character of the acoustic field inside a wind tunnel, can be used. The model 
consists of a known harmonic acoustic source of finite dimension inside an infinite 
duct of constant cross-sectional area. The duct is assumed to contain a uniform 
subsonic flow. The wall boundary will be represented by a simple impedance 
condition. 

This simpler model assumes that several physical effects are negligible. First, 
no wind tunnel is an infinitely long uniform duct; all contain drive fans and 
vanes, changes in area and shape, and are closed loops or open loops. All of 
these factors will affect the transmission of acoustic energy in the wind tunnel. 
Often, these changes are far from the test section which contains the noise 
source and the acoustic field of interest; therefore, it can be assumed that the 
changes have a negligible effect on the sound field in the test section. Vanes 
will have a small amount of acoustic backscatter. The drive fan will be far 
from the test section so its presence as an extra noise source and its influence 
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in scattering sound should be small. Most changes in the cross section of a 
wind tunnel are gradual so they will not produce a strong acoustic reflection. 
The topological departure from an infinite duct is manifested by changes in the 
geometry far from the test section. The fundamental acoustic mode of the wind 
tunnel will be a very low frequency since the wavelength will be of the order 
of the length of the tunnel in a closed circuit tunnel, or twice the length of 
the tunnel in an open tunnel. High-frequency structual modes of the tunnel, 
where the acoustic source frequencies occur, should not strongly couple into the 
problem. 

Next, the boundary conditions at the wall are not exact It is assumed 
that the influence of the boundary layer is small and can be incorporated into 
a point impedance condition. This local impedance condition also precludes 
any coupling between the acoustic field and structural modes of the wind 
tunnel. 

Last are the assumptions involving the source This model assumes that the 
source radiation is known and that the amplitude is low enough for the linearized 
acoustic equation to apply For convenience, it is also assumed that the source 
can be decomposed into its separate frequency components, implying that the 
source is both periodic and that the linear equation holds Finally, it is assumed 
that the presence of the walls does not affect the acoustic radiation from the 
source. This is reasonable, because the acoustic source is not ty^pically close to 
any wall; as a result the pressure and velocity induced by the presence of the 
wall are small compared to the pressure and velocity at the source produced 
by the noise sources on a helicopter These assumptions are consistent with 
modeling the source as a distribution of rotating acoustic monopoles, dipoles, and 
quadrupoles 


MATHEMATICAL FORMULATION OF PROBLEM 


An inhomogeneous, hyperbolic partial differential equation with mixed 
boundary conditions on the walls models the general problem of an acoustic 
source inside an infinite duct with uniform subsonic flow No analytic solution 
is known for this problem with arbitrary sources, arbitrary duct shape, and 
arbitrary wall impedance In this section, the problem will be formulated as a 
differential equation with boundary conditions and transformed into an integral 



equation. 


The acoustic problem is modeled with linearized potential flow. Acoustic 
pressure P will be the variable used. A complex pressure is used to keep track of 
the phase information The physical acoustic pressure is represented by the real 
part of P. 

The inhomogeneous convected wave equation with appropriate boundary con- 
ditions describes this acoustic flow. The wave equation comes from consider- 
ing a small perturbation to the Euler equations describing flow in an inviscid, 
homogeneous fluid which does not conduct heat. A combination of the equations 
describing conservation of mass and momentum will produce the wave equation, 
as in reference 7 For a uniform mean flow &oo with acoustic source term F the 
governing equation is 

V 2 P(?, t) - Uj- t + &oo • V) 2 P(r, () = F(r, t) (1) 

The boundary condition on the walls of the duct relates the acoustic pressure to 
the normal gradient of pressure: 


aP — /?VP • n = 0 (2) 

The coefficient combination ft /a. has the dimensions of pressure/velocity and is 
equivalent to specifying the impedance without flow The normal to the duct 
wall is n (positive inward) Also, a radiation condition must be enforced in the 
duct far away from all sources The radiation condition used will be explained 
later when the integral relation is developed for the acoustic potential in a 
duct. 

The coordinate system, as shown in figure 1, is oriented such that the only 
mean flow component is in the positive x-direction. No loss of generality results, 
because the flow is assumed to be uniform in a duct of constant cross section. The 
governing equation becomes 


V 2 P(f, t) - 1(|- + Uoo ■ |-) 2 P(r, t) = F(r, t) 
c 2 at ox 


( 3 ) 


5 



T 

n y 

* 


Uoo, M 




Figure 1. - Fixed reference frame 


The acoustic sources considered and, therefore, the resulting acoustic field, 
will be periodic. A Fourier series decomposition will be used so each frequency 
component can be considered separately. The acoustic pressure and source will 
now have the forms 

P(G *) = ]£ a(w)e _tw *p(r) (4) 

u 


and 

F(?,t) = E i ’( w ) e_ ”'/( 7 > (5) 

6 0 

With this transformation there is a change to the variable k = u/c. The equation 
becomes 

V 2 p(?) - M 2 + 2iMk + k 2 p (?) = /(?) (6) 

The free-field Green’s function for the convected Helmholtz equation will 
be used. The Green’s function, Gh(r;r s ), describes the acoustic response at 
the location ? to a harmonic source of frequency u at the location f s in 
unbounded space For a uniform mean flow in the positive x-direction, Gh(r,r s ) 
satisfies the following partial differential equation with a radiation condition at 
infinity: 
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V 2 G h (?;?,)-M 2 ^i^+2«m^^l+*: 2 G h (?;? s ) = « 3 (?-? s ) (7) 


The solution (derived in appendix A) is 


Gh (r;r s ) 


(1— M 2 ) 


4ttR 


(8) 


with 

R = \j(x - x. ) 2 + (1 - M2)[(J, - y, )2 + (Z - 2 t )2] (9) 

In unbounded space, this Green’s function gives the solution to equation (6) by 
evaluating the integral 


OO OO 00 

P/(r)= / / / G h( r;r s )f(? s )d? s (10) 

— 00 — 00 —oo 

The subscript f on the pressure coefficient will be used to designate this pressure 
as the pressure in unbounded or free space. Next, an integral relation is formed 
for p(r) on and inside the duct by use of Green’s second formula The relations 
are 


- IT = / / (-p(?)[ VG h(r;r s ) «] + ?*)[Vp(r) • n])dS + (11a) 


til 


G h (r,r s )f(r)dV 


on the duct surface S, and 

V{r) = J J (-p(?)[VG h (r;r s ) n] + G h (f, f,)[Vp(r) • n])dS + (116) 


/// 


Gh(f;r s )/(7)dV' 
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in the duct volume V. The derivation of this relation is given in appendix B. 
The corresponding equations for the two-dimensional case are given in appendix 
C. 


On the surface of the duct S, the integral relation and the boundary condition 
combine to give 


pQ 

2 


// 


a 


P(?)|- V G h (r ; r s ) • n + Gh ( r ; r s )- j 


dS + 


( 12 ) 


/// G h (?;r,)/(r)dy 

when (3^0. The restriction that is applicable to any solid or porous surface 
since ft = 0 implies a vacuum. This is an exact equation for the acoustic 
pressure p(r) on the surface of the duct. All quantities except p(r) are known. 
In general, this equation cannot be solved analytically; however, it can be solved 
numerically. From the solution of p(r) on the surface of the duct, the value of 
p(r) anywhere inside the duct, and thus the acoustic pressure in the field, can be 
found from 


P(r) 


= // p(?) ~ 


a 


VGh(r;r s ) n + G h (r;r s )- 


d£ + 


(13) 


III Gh ^;r s )f(r)dV 

In equations (12) and (13), the term / / /Gh(r,r s )/(r)dY is equal to p^(r), 

the acoustic pressure which would exist with the same source, f(r), in uniform 
flow and no walls. This term can be calculated in many ways. One method is 
described in the following paragraph. 

It is convenient to calculate the acoustic pressure field in the frequency domain 
since the reverberant field will be calculated in the frequency domain. The follow- 
ing formulas give the acoustic pressure spectrum for the rotor mth blade passage 
harmonic, u = mNO. The formulas were taken from reference 8; the observer 
is in motion to give the wind-tunnel reference frame. Subcript h refers to the 
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hub: 


P f(x,V,z) = 


mN 2 fi 


47rR[l — Mcos(<5 r )] 


♦mNfl[R — M(s — xh )] 

e (i-m2) c 


n= ^oo J o 'Vc[l-Mcos((5 r )]; 


mftyo(mN — n)i(mN— n)V n -f- ^®l( m N-n) dr 


The coefficient of the thickness term is 


/*£*« *nf 

t n= L f (^ r ) e r d £ 

''{le 


The coefficient of the loading term is 


tn€ 

U = / /(f,r)e r d£ 

J tie 


The other variables are. 


(14) 


(15) 


(16) 


R 


)/(*— z/i) 2 + (1 — M 2 )[(2/— y^) 2 +( 2 — ^) 2 ] 


(17a) 


and 


1 Mcos(tf f )= ^ M ^ R 
v ' R — M(z — x s ) 

(176) 

v. St*:! 

(17c) 


The acoustic pressure in the wind tunnel (duct) is found by solving the 
equations in several steps. First, the free-field pressure p^ is found on the inside 

surface of the duct. Evaluation of equation (14) is one example of how this 
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term may be found. Substituting p f = / // Gh(r;r s )/(r)dV into equation (12) 
yields an integral equation for the acoustic pressure on the inside duct surface. 
Solution of the integral equation will yield the acoustic pressure on the duct 
surface. Then, the acoustic pressure anywhere inside the duct may be found by 
substituting the surface pressure into the integral in equation (13) and evaluating 
the integrals. 


10 



APPENDIX A 


DERIVATION OF THE FREE-FIELD GREEN’S FUNCTION FOR THE 
SUBSONIC CONVECTED ACOUSTIC EQUATION 


The following partial differential equation (Al) and boundary condition (A2) 
for G(r, t; r$,t s ) describe the effect of an acoustic source at the location 7 S with 
time dependence f(t s ) on a receiver at the location 7 and time t in unbounded 
space A uniform subsonic flow with Mach number M ip the positive x-direction 
is modeled. 


( v2 - 


M' 


d 2 2M d 2 


dx 2 


c dxdt 


\_cP_ 
c 2 dfi 


) 


G(?,t,r s ,t s ) = f(t s )6z(?-? s ) (Al) 


The boundary condition at infinity requires a radiation condition at infinity 
of 


lim 

?-+oo 


(7-7 s )(7 VG *(?;?,) + 


dGh^7s) ] 

c dt } 


= 0 


(A2) 


This ensures that waves are outgoing at large distances and that contributions to 
surface integrals at infinity are zero 


In the original coordinates (fig 1) the observer and source, 7 and 7 S , 
are fixed in space and the fluid is moving in the positive x-direction with 
Mach number M. A solution can be found with a Galilean transformation 
to the reference frame where the fluid has no time-averaged velocity (M = 
0) and the source and observer are moving uniformly in the negative x- 
direction with Mach number M. The transformation shown in figure 2 is as 
follows: 


r = ( x,y,z ) 

(A 3a) 

7 = r — Met 

(A 3b) 

x = x — Mcf 

(A 3c) 

y = y 

(A 3d) 

z = z 

(A Se) 
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t = t 


(AS/) 



Figure 1. - Fixed reference frame. 



r s (t > 0) r s (t - 0) 


n 



Figure 2. - Moving reference frame. 
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Equation, Al, transforms to 


~2 

V 


c 2 at 2 


G (T/t;7 s ,i s ) = f(i s )6z{? - ? s ) 


(A4) 


The transformed equation has a solution, which is an integration over time and 
space, 


G (7,t;f s ,i s ) = 


(AS) 


00 oo oo oo 

1 1 1 1 

— oo —oo —oo — 00 


f{ts)h{^-7s)8[ t s 


(■-*)) 


47t|? — 3| 


dt s d3 


A spatial integration with the three-dimensional delta function leads to a time- 
integration, 


00 

G {7,t,? s ,is)= J 

— 00 


fO. )s 


(M‘~ 


I r—r. 


4tt| 7 — ? s 


-At* 


(A6) 


A time-integration with the delta function yields 


G(?,'t;? t ,k) 




\ 

— /(*«) 




M 

J 


(Al) 


at t$ — t — 


I*-*. I 


The moving source and observer have the following forms, a constant plus a 
component linearly proportional to time: 


?g 7 go 


— Mct s 


(A8a) 
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Xg Xgo Mcfg 

y s = y S o 

Zg ZgQ 

f = ? 0 — Me t 



y — Vo 
z = z 0 


I?— ?s| = \] (£o — Me t — iso + Mc^s ) 2 + (y 0 — y so ) 2 + (20 — ^so ) 2 
A source-time derivative of the distance is 


_d_ 

dt s 


|?-?s| = 


Mc(io — Met — x so + Mcf s ) 


\?-?s\ 


or 


_d_ 

dtg 


\r-f a \ 


Mc(x 0 — Zso) J^j 2 c 

\?-?s\ 


Evaluating the derivative gives 


G(f, t;?s, ~t s ) 



M, | v 

C ) 


47r[|f-? s |(l-M 2 ) + M(5 0 



In the moving coordinate system, 


G(?,t;?s,~ts) 



M.l 

c / 


4tt[( 1 - M 2 )\7 - ? s | + M(x - x s + Met - Mcf s )] 


(A8b) 
(A8c) 
(A8d) 
(A 8e) 

(A8f) 
(ASg r) 
(A8h) 

(A St) 
(A9a) 

(A9b) 


(A10) 


(All) 
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Now the original coordinates from the original reference frame can be sub- 
stituted into the formula for G(r, t, r s , t s ). 


G(f, t;r 8 ,t t ) 



|F— $lct— f e +Klct s | 
c 


) 


47 t [(1 — M 2 )\? — filet — ? s + filct s | + M(x — x s )] 


(A 12) 


c(t — t s ) = \r — Met — r s + Mct s | (A13) 

This equation needs to be solved for the source time so that the formula for 
G(r, t; r s , t s ) can be evaluated: 

c 2 (t — tg) 2 = (x — x s — Mcf-|- Mct s ) 2 + (y ~ Vs) 2 +(2 — z s ) 2 (A14) 


= {x— x s f +2(a: — x f )Mc(f s — t) + M 2 c 2 (£ — i s f +{y— y s ) 2 + (2 — z s ) 2 
Collecting terms yields a quadratic equation in c(£ — t s ): 

(1 — M 2 )c 2 (£ — t s ) 2 + 2M(z — x s )c (t — t s ) — (A15) 


(a: - a: s ) 2 - (y - y s f - {z - z s f = 0 


with solutions 


c (t- t s ) = 


— M(z — x s ) 
1 - M 2 


(A 16) 


y/M 2 (a:— j s ) 2 +(1 — M 2 )[(a:- z s ) 2 + {y — y s ) 2 + (z — ^s) 2 ] 
± 1 - M 2 


The positive root is used to ensure causality 


R = 


— M(z — z s ) , 

1 - M 2 + 


(All) 


y/M 2 (a:— x s ) 2 +(1 — M 2 )[(a: — z s ) 2 + (y — t/ s ) 2 + (z — £ s ) 2 ] 

1 — M 2 
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Substituting the solved value of t s gives 


G(r, t; r s ,t s ) 


4tt[( 1 — M 2 )R + M(a: — a: s )] 


(A 18 a) 


or 


G (r,t;r s ,t s )= (AlSb) 

t( + ~ M(s— s 8 )+ ^/M 2 (*— T e ) 2 +(1— M 2 )[(*—*» ) 2 +(y— V» ) 2 +(«— ) 2 1 A 

Jy c(i — m 2 ) J 

A time-dependent harmonic source is substituted into the equation. Let 
f(t ) = e~ wt ,k = then, the partial differential equation becomes 

(V 2 — M 2 — (- + k 2 )G h (r,r s ) = fa(r — r s ) (A19) 

dx 2 dx 

and the solution becomes 


G fc (r;? f ) = 


jf +\Z(j— ag) 2 +( 1 ~ M 2 )[(v— ys) 2 4-(z— ^s) 2 ] — M(a— acs) ^ 

_ e v y 

4 tt\/(^ — £s ) 2 +(1 — M 2 )[(t/ — y~ s ) 2 + (2 — ~ z 7 F ) 


(A20) 


This solution gives waves that propagate away from the source and 
decay with distance. Except in the case of M = 0, the surfaces of 
constant phase and constant amplitude do not coincide. The function 
G/i(r;r s ) is very well behaved far from the source and has the following limit 
properties: 

lim Gh(r;r s ) = 0 (A21) 

f—*oo 

dGh (?;? s ) n „ - , ^ 

lim - = 0, ^ = ax + by + cz 

?-* oo ds 
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lim 

x-*oo 


lim Gh(r;r s ) = 0 

x — >00 


dG h {7;7 s ) 

ds 


s = ax + by + cz 
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APPENDIX B 


DERIVATION OF GREEN’S FORMULA FOR THE SUBSONIC CONVECTED 
HELMHOLTZ EQUATION IN A UNIFORM DUCT 


An integral equation can be constructed which describes the harmonic 
acoustic pressure coefficient p inside an infinite duct with a uniform sub- 
sonic flow by utilizing the Green’s function for the acoustic pressure and 
Green’s second formula. The subsonic flow in the positive x-direction has 
Mach number M. The equation for the Green’s function in subsonic flow 
is 

(V 2 - M 2 + 2tkM^- + k 2 ) G*(r; f, ) = S 3 (r - r , ) (fll) 

dx 2 dx 

And the equation for the acoustic pressure from an arbitrary source is 

(V 2 — M 2 + fc 2 )p(r) = /(r) (B2) 

dx 2 ax 


By restricting the formulation to cases in which a/ [3^0 and a/^oo, it can be 
shown from the differential equation formulation of the problem that the acoustic 
pressure in the duct will decay with distance away from the source Thus, the 
far-field pressure has the following form: 


lim p(r) = 0 (BSa) 

f— *oo 


dp(?) ^ 

lim — — 1 = 0, s = ax + by + cz 
?— >oo as 


Green’s second formula is 



AGh(r; 7 s) 

dn 


H~ G/ l (?*;? s ) 


<9p(?) 

dn 


dS = 


(BSb) 


(B 4) 


18 



///[ 


p(?)V 2 G A (?;? s )- G*(?;? s )V 2 p(?) 


dV 


with the surface integral evaluated over the inside surface of the duct and the 
volume integral evaluated over the interior of the duct. The normal to the surface 
with positive direction inward is n. Substituting V 2 G/ l (f;? s ) and V 2 p(?) into 
equation (B4) gives 


//( 




dn 


dn 


dS = 


(B 5) 


III™ 


r 2 d 2 G h (?;? s ) 


■9G h (f;? s ) . 2i 


+ 2iA:M— - k* G fc (7; r s ) + <5 3 (? - r s ) 


dV- 


///cm?;?,) 


M‘ 


«9 2 p(y > ) 

da; 2 


r W) 


+ v - — fc 2 p(r) + /(?) 

da; 


dV 


Canceling terms gives 


//[• 


■m5#l +wlll W) 


dn 


dn 


dS = 


(B 6) 


III 

III 


A (r,/«) a2p(7) 


da; 2 


da; 2 


dV + 




da; 


da; 


dV + 


III W)W-?s)-G h (7;? s )f(7)]dV 


The terms in the first integral on the right-hand side of equation B6 
can be evaluated analytically by integrating by parts for a; = — oo to 
-(-oo 


/// 




dV = 


(B7) 


19 



// 
-III 


M»f P (r)- 9G * (? * ) 


dx 


~Gh(?; 7 s) 


w) 

dx 


-*\ 1*=00 


dA 


■ix =— 00 




2 r dm dGm f_s ) <9Gfr(7; ) flp(7) 

dx dx dx dx 


dV = 0 


with the integration over the cross-sectional area of the duct in the surface 
integral. Both integrals are identically zero. The terms in the area in- 
tegral, 


r dGhir'As) r, r - ^p(r) A 
p (r) — G h (r;r s )-£A = 0 


(£ 8 ) 


at x = ±oo, since as x-+oo, G/*(r;r s ) = 0, = 0, p(r) is finite, and 

is finite. 

The second volume integral in equation B6 is also identically zero. Again 
integration by parts produces terms that are all zero: 


m\ 




dx 


dx 


dV = 


IB 9) 


III 




dx 


dA- 


— 00 


/ / jH Un a ^^-G h (f,r,f^ 


-oo 


dx 2 


dx 2 


dzdA = 


Il\^f«>*-TI™A"' 


dA = 0 


Both the functions G/ l (r; r s ) and p(?) are well behaved enough so that at 

x = ±oo, = 0, = 0, / p(r)da; is finite, and / Gfc(?;? s )da: is finite. 

Applying these constraints, an integral relation is found for the acoustic pressure 
p(r) on the surface, 


P(r) [ f\ ^ dG h(r;r s ) , ^ M?) 


= //h 


P(?)- 


dn 


+ Gh(7*;? s )- 


dn 


dS + 


(£10a) 
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Ill G ^;r s )f(r)dV 
And in the duct volume the equation is 

p(7) _ / / [.^ + G*(f;r.)^]dS + (Bm) 
-/ / / G*(?;?,)/(?)dV 
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APPENDIX C 


INTEGRAL EQUATION FOR TWO-DIMENSIONAL DUCTS 


The three-dimensional case can easily be reduced to a two-dimensional 
problem by assuming that all quantities are independent of z. The walls 
become two parallel plates at y = t/j and y = y u and the sources be- 
come infinite line sources in the z-direction. The resulting acoustic field is 
also uniform in the z-direction. Putting these restrictions into the three- 
dimensional equation simplifies the equation so that Gh(r;r s ) is the only 
quantity that depends on z. The integrals in the surface and volume equa- 
tions take on the following form with the restrictions that and 

oi/f3^oc: 

OO 00 

/ J V{x,y\,^-^ + ^^^G h {r;r s )dzdx+ (Cl) 

— 00 — OO 


00 OO 

J v{x,y u ,z)(^^ + ^^^Gh(r;r s )dzdx + 

— 00 —00 

oo oo 

/ / / f(x,y, z)G h (7,? s )dzdxdy = 

S/j —oo — oo 

00 00 

/ p(z> 2/l)(- ^ / G/ l (r;r s )d^do: + 

— OO — 00 




G/i(r, ? s )d2da; + 
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If)l 00 00 

II /(*, y) i Gh{r,r s )dzdxdy 

S/j —oo —oo 

The integration in the z-direction can be done analytically (ref. 9) for all of the 
integrals since 

OO 

G 2 d{x,y;x s ,y s )= j Gfc(r;r s )d2 = (C2) 

— OO 


/ 


— 00 


.i.( / (*— X S ) 2 +(l— M 2 )[(y— Vs) 2 +(«— 2t) 2 ]- M(*— x s 

\ (1 M 2 ) J 


4 tt\/(x — z s ) 2 +(1 — M 2 )[(y — i/ s ) 2 + (2 — 2 S ) 2 ] 


d z = 


y/1 - M 2 e- tkM ( x ~ x °) . ( ky/(x - z s ) 2 + (1 - M )(</ - ?/ s ) 2 


4 ki 


Hjf 


1 - M 


) 


This provides a two-dimensional Green’s function for the convected Helmholtz 
equation which will be called G 2 d{x,y,x s ,ys) The function Hq is the Hankel 
function of the first kind of order 0. The surface integral equation is 
now 


p(z, y) 
2 


(C 3) 


/ 


— 00 




y\,x s ,y s )dx + 


S 


— OO 


( ft fytkc i 

— + ~p~jG 2 d(x, y u ; Xs , y s )da: + 


S^u oo 

/ / 


J/J — oo 


f(x, y)G 2 d(x, y; x s , y s )dzd?/ 
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This is now the equation to be solved for p {x,y). All quantities are known 
except p(a;, y). As before, the last integral represents the complex acoustic 
pressure that would exist in free space without the duct from the source 
f(x,y ) and can be replaced by Pf(x,y) if the free acoustic field is known. 
After p(x, y) is solved for on the surfaces y u and y\ the acoustic pressure 
and thus all of the acoustic quantities, can be evaluated inside the duct 
with 

p(x, y) = (C4) 

00 

J p(z, 2/i)(- ^ + ^j^jG 2 d(x, y\; x s , y s )dx + 

— OO 


/ 


— OO 




x, y VL ;x s ,y s )dxdy -\- 


y u oo 

/ / 


yj — oo 


f(x, y)G 2 d(x, y, x s , y s )dxdy 
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